A SPARSE FLAT EXTENSION THEOREM FOR MOMENT 

MATRICES 



MONIQUE LAURENT AND BERNARD MOURRAIN 



Abstract. In this note we prove a generalization of the flat extension theorem 
of Curto and Fialkow |4] for truncated moment matrices. It apphes to moment 
matrices indexed by an arbitrary set of monomials and its border, assuming 
that this set is connected to 1. When formulated in a basis- free setting, this 
gives an equivalent result for truncated Hankel operators. 



1. Introduction 

Throughout this note, K denotes a field, IK[x] = IK[xi, . . . , x„] is the ring of 
multivariate polynomials in n variables x — (xi, . . . , a;„) with coefficients in K, 
Mn — {x" :— ■ ■ ■ a;"" | a G N"} is the set of monomials in the variables x, and 
■Mn,t (resp., IK[x]t) is the set of monomials (resp., of polynomials) of degree at most 
t. The dual basis of Mn in the dual space K[x]* is denoted as P„ = {d'' | /3 G N"}. 
The natural action of K[x] on K[x]* is denoted by 

{p, A) G ]K[x] X IK[x]* ^p-Ae K[x]* 

where {p ■ A){q) A{pq) for q G IK[x]. 

1.1. The moment problem. In this section, we consider IK = R. The moment 
problem (see e.g. [B Ej) deals with the characterization of the sequences of mo- 
ments of measures. Given a probability measure fi on R", its moment of order 
a = x" G Ain is the quantity Jx°'fi(dx). The moment problem concerns the 
characterization of the sequences y — {ya)aeMn that are the sequences of moments 
of some nonnegative measure /i, in which case one says that /i is a representing 
measure for y, with yi = 1 if /i is a probability measure. Let A G R[x]* denote 
the linear form on ]R[x] associated to the sequence y, defined by A{p) — J^aPaUa 
for any polynomial p = X^aeAi P"-'^ ^ M[x]. Then, y has a representing measure 
fi precisely when A is given by A{p) — J p{x)fi{dx) for all p G R[x]. A well known 
necessary condition for the existence of a representing measure is the positivity of 
A, i.e. A(p^) > for all p G M[x], which is equivalent to requiring that the matrix 
M{y) := {yab)a.beM„ be positive semidefinite. As is well known this necessary con- 
dition is also sufficient in the univariate case {n = 1) (Hamburger's theorem), but it 
is not sufficient in the multivariate case (n > 2). However, positivity is sufficient for 
the existence of a representing measure under some additional assumptions. This 
is the case, for instance, when the sequence y is bounded [2] or, more generally, 
exponentially bounded [3]. The next result of Curto and Fialkow (4J shows that 
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this is also the case when the matrix M{y) has finite rank (cf. also [IH [15] for a 
short proof). 

Theorem 1.1. [4] IJ M[y) is positive semidefinite and the rank of M(y) is finite, 
then y has a (unique) representing measure (which is finitely atomic with rank M{y) 
atoms). 

In the univariate case n — \, a. matrix of the form M{y) is a Hankel matrix. In 
the multivariate case, M{y) is known as a generalized Hankel matrix (see [17]) or 
moment matrix (see |15j). One can also define truncated vaomcni matrices: A matrix 
M indexed by a subset C C is said to be a moment matrix if Ma^b — Ma'y 
for all a,b,a',b' £ C with ab ~ a'b'. Thus its entries are given by a sequence 
y = (2/c)cgC-Ci where C ■ C := {ab | a, 6 G C}, and we can write M = Mc{y). 
When C = A4n.t, we also write M = Mt{y), where the entries of y are indexed 
by A^n,2t- Such matrices arise naturally in the context of the truncated moment 
problem, which asks for the existence of a representing measure for a truncated 
sequence indexed by a subset of monomials. A solution to the truncated moment 
problem would in fact imply a solution to the moment problem. Indeed, Stochel 
[iS] shows that a sequence y = (ya)aex„ has a representing measure if and only if 
the truncated sequence (ya)aeA<„ t has a representing measure for all t gN. 

1.2. The flat extension theorem of Curto and Fialkow. Curto and Fialkow 
studied intensively the truncated moment problem (cf. e.g. [U [Sj [6] and further 
references therein). In particular, they observed that the notion of flat extension of 
matrices plays a central role in this problem. Given matrices Mq and Afg indexed, 
respectively, by C and B Q C, Mq is said to be a fiat extension of Afg if Afg coincides 
with the principal submatrix of Mq indexed by B and rank Mq — rank M^. Curto 
and Fialkow [4J show the following result for truncated moment matrices. 

Theorem 1.2 (The flat extension theorem [4]). For a sequence y = {ya)aeMn 2tr 
if Mt{y) is a flat extension of Mt-i{y), then there exists a (unique) sequence y = 
{ya)aeM„ fof which M{y) is a flat extension of Mt{y). 

The flat extension theorem combined with Theorem 11.11 directly implies the 
following sufficient condition for existence of a representing measure. 

Corollary 1.3. For a sequence y — {ya)aeM„ 2t ^ ¥ Mt{y) is positive semidefinite 
and Mt{y) is a fiat extension of Mt-i{y) , then y has a representing measure. 

Curto and Fialkow [5] show moreover that the flat extension condition is in some 
sense necessary and sufficient for the existence of a representing measure. More 
precisely, they show that a sequence y = {ya)a£Mn 2* has a representing measure if 
and only if it can be extended to a sequence y' = iy'a)aeM„ 2t+2fc+2 (fo^' some fc > 0) 
for which Mt+k+i{y') is a flat extension of Mt+k{y')- 

The proof of Theorem [TT2] relies on a "truncated ideal like" property of the kernel 
of flat moment matrices (see (|2.2p below). This permits to set up a linear system 
of equations in order to construct the flat extension Mt+i{y) of Mt{y) (and then 
iteratively the infinite flat extension M{y)). See also [T3] for an exposition of this 
proof. Schweighofer [18] proposes an alternative proof which is less technical and 
relies on properties of Grobner bases. We propose in this note another simple alter- 
native proof, which applies more generally to truncated moment matrices indexed 
by (suitable) general monomial sets (see Theorem II. 4|) . 
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1.3. A generalized flat extension theorem. We need some definitions to state 
our extension of Tlicorcm ll.2l For C C A4n, 



n 



:= CU (J XiC — {to, xim, . . . , Xnm \ m £ C} and dC := \ C 

i=l 

are called, respectively, the closure and the border of C. The set C C Mn is said 
to be connected to 1 if 1 € C and every monomial to S C \ {1} can be written as 
m — Xi-^ ■ ■ ■ Xi^ with , Xij^Xij , . . . , • ■ ■ Xi^ G C. For instance, C is connected to 
1 if C is closed under taking divisions. For example, {1, X2, xiX2\ is connected to 1 
but {1, XiX2\ is not. We now state our main result. 

Theorem 1.4. Consider a sequence y = (2/a)aGC+ c+; where C C M.^ is finite and 
connected to 1. If M^-^-iy) is a flat extension of Mc{y), then there exists a (unique) 
sequence y = {y)aeM„ for which M[y) is a flat extension of M^+iij). 

The proof is delayed till Section [2l Note that Theorem 11.21 follows directly from 
Theorem 11.41 applied to the case C — M.n,t-i- Thus our result can be seen as a 
sparse version of Theorem 1 1.21 which applies to a more general monomial set C, not 
necessarily the full set of monomials up to a given degree. We now give an example 
showing that the assumption that C is connected to 1 cannot be omitted. 

Example. For n — 1, consider the set C — {1, x"^}, which is not connected to 1, with 
dC = {x^x"^}. Consider the sequence y S M''^ ''^ defined by yi = y-x = Vx^ = 1, 
Ux^ — Dx^ = J/x's ~ o, and yx^ — yx'' = Vx^ — where a, b are scalars with b ^ a^. 
Then, rank Mc+{y) — rank Mc(y) — 2. If there is a flat extension M{y) of Mc+{y), 
then its principal submatrix indexed by C"*" U {x^} has the form: 



Mc+ulx-^jiy) 



However, 1 — x G ker Mc+{y) implies x — x'^ £ ker Mc+u{x^}{y) (see p.2p ) and thus 
1 = a = b, contradicting our choice b ^ a^. Hence no flat extension exists. 

1.4. Basis- free reformulation. Here we reformulate our result in a basis-free 
setting. Moment matrices correspond indeed to choosing the monomial basis Mn 
in the polynomial ring K[x] and its dual basis Vn in the dual space K[x]*. Given 
A e K[x]*, the operator 

Ha ■■ K[x] K[x]* 
p p ■ A 

is known as a Hankel operator. Its matrix with respect to the bases and 
Vn is precisely the moment matrix {h.{x"'^^))a,i3eM„ = ^Hv) of the sequence 
y = {A{a))a(zMr,- The kernel of Ha, 

ker Ha = {p£ ]K[x] | A{pq) = G K[x]}, 

is an ideal in K[x]. Moreover, when K = M and A is positive, i.e. when A(p^) > for 
all p € K[x], ker Ha is a real radical ideal [2]. Theorem 11.11 means that A G M[x]* 
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is positive with rank Hf^ < cx) if and only if there exists a nonnegative finite atomic 
measure /i for which A(p) = j p{x)iJL{dx) for aU p G R[x]. 

Truncated Hankel operators can be analogously defined. Given C C A1„ and 
A G (Span(C"'" ■ C+))*, the corresponding Hankel operator is 

Hi* : Span(C+) Span(C+)* 
p ^ p ■ A 

and its restriction to Span(C) is : Span(C) Span(C)*. We have the following 
mappings: 

Span(C) ui Span(C) 0-2 Span(C+) 



1.1 



ker Hi ker H^^ n Span(C) ker 



where cti is onto and (72 is one-to-one, so that 

, ,. Span(C) Span(C) Span(C"'") 

1.2) dim — — V < dim < dim — — 

^ ^ ker Hi " ker Hl^ n Span(C) " ker 

Thus, rank Hi* — rank Hi (in which case we also say that Hi* is a flat extension 
of Hi) if and only if equality holds throughout in (|1.2p . i.e. both ai and (T2 in (|l.ip 
are isomorphisms or, equivalently, if 

Span(C+) = Span(C) + keri7^^ and ker Hi = ker Hi* n Span(C). 
Theorem 11.41 can be reformulated as follows. 

Theorem 1.5. Let A G (Span(C^ -C^))*, where C C A4„ is finite and connected to 
1, and assume that rank Hi = rank Hi. Then there exists (a unique) A G R[x]* 
for which Hf^ is a fiat extension of Hi* , i.e. A coincides with A on Span(C"'' ■ C^) 
and rank Hj^ — rank Hi* . 

1.5. Border bases and commuting multiplication operators. We recall here 
a result of [TB] about border bases of polynomial ideals that we exploit to prove our 
flat extension theorem. Let B := {&i, . . . , b^} be a finite set of monomials. Assume 
that, for each border monomial Xibj G dB, we are given a polynomial of the form 

N 

gi^i) ■= x,bj - J2 at'^^h where a|;^^ G K. 

h=l 

The set 

(1.3) F := | i = 1, . . . , ?i, j = 1, . . . , TV with x,bj G dB} 

is known as a border prebasis [S] or a rewriting family for B [12] ■ When the set B 
contains the constant monomial 1, one can easily verify that i3 is a generating set 
for the quotient space IK[x]/(i^), where (F) is the ideal generated by the set F. 
When B is connected to 1, Theorem II .61 below characterizes the case when i3 is a 
basis of K[x]/(F), in which case F is said to be a border basis of the ideal (F). For 
this, for each i = 1, . . . , n, consider the linear operator: 

Xi ■ Span(;B) Span(B) 

extended to Span(;B) by hnearity. When ;B is a basis of K[x]/{F), Xi corresponds 
to the "multiplication operator by Xi" from K[x]/(F) to IK[x]/(F) and thus the 
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operators commute pairwise. The next result of |16j shows that the 

converse implication holds when B is connected to 1; this was also proved later in 
[5| when B is closed under taking divisions. 

Theorem 1.6. 16J Let B C A4„ be a finite set of monomials which is connected 
to 1, let F be a rewriting family for B as in (|1.3p . and let Xi, . . . ,Xn be defined 
as in (|1.4p . The set B is a basis of the quotient space K[x]/(_F) if and only if the 
operators xi, • ■ • , Xn commute pairwise. 

The proof of our sparse flat extansion theorem is an adaptation of this result to 
kernels of Hankel operators, where we omit the assumption that B is connected to 
1. 

1.6. Contents of the paper. Section [2] contains the proof of our generalized flat 
extension theorem and we mention some applications in Section[3l In particular, we 
observe that Theorem II .21 is an 'easy' instance of our flat extension theorem (since 
one can prove existence of a basis connected to 1). We also point out the relevance 
of the flat extension theorem to polynomial optimization and to the problem of 
computing real roots to systems of polynomial equations. 

2. Proof of the flat extension theorem 

We give here the proof of Theorem 11.51 (equivalently, of Theorem 11.41) . We 
will often use the following simple observations, which follow directly from the 
assumption that rank ^ rank H^: For all p G Span(C^), 

(2.1) p e ker H^^ 4^ A{ap) = Va e C+ A{ap) = Va e C, 

(2.2) pekeiH^^ and x^p e Span(C+) Xip <E ker H^^. 

Our objective is to construct a linear form A G ]K[x]* whose Hankel operator if^ is 
a flat extension of i/^^ . 

Let B C C for which rank = rank = Note that we can assume that 
1 £ B. Indeed, if no such B exists containing 1, then A(p) = Vp £ Span(C"'") and 
one can easily verify that this implies that A is identically zero, in which case the 
theorem trivially holds. 

From the assumption: rank H^^ = rank = \B\, we have the direct sum 
decomposition: Span(C"^) Span(S) © ker , and thus 

(2.3) Vp e Span(C+) 3! n{p) £ Span(S) such that f{p) :=p-7r(p) £ ker H^'' . 
Then the set 

F := {f{m) ~ m ~ 7r(m) | m £ dB} 
is a rewriting family for B and, for i = 1, . . . , n, the linear operator Xi in (11.4p maps 
p £ Span(S) to Xiip) = T^i^ip) £ Span(;B). We show that xij---,Xn commute 
pairwise. Set K := ker . 

Lemma 2.1. Xi ° Xj = Xj ° Xi- 

Proof. Let m £ B. Write n{xim) := J2bei3K^ (K *= We have: 

Xj o Xiim) = XjiJ2beB ^h^) = EbGB KXjib) = J2beB Ki^jb - f{xjb)) 
= XjiJ2beB ^6^) - T.beB Kfi^jb) = Xj{x,m - fix.m)) - XibeB Kfi^j'^)- 
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Therefore, 

' " ' beB 

Pi , ■' 

P2 

We show that pi G K. Indeed, Va G C, A(api) — A{axif{xjm) — axj(xim)) — 
since axi,axj € and /(x^m), f(xjm) e by (|2.ip . this shows that pi G if. 
As p2 G ^ too, this imphcs p E K and thus p = 0, because p € Span(S). □ 

Our objective now is to show that ;B is a basis of K[x]/(F) and that, if tt denotes 
the projection from K[x] onto Span(;B) along (F), then the operator A defined 
by A{p) = A{7t{p)) for p G K[x], defines the desired flat extension of A. Note 
that when B is connected to 1, Theorem 11.61 implies directly that S is a basis of 
K[x]/(F). As we do not assume B connected to 1, we cannot apply Theorem 11.61 
but our arguments below are inspired from its proof. In particular, we construct 
the projection tt via the mapping (p from (|2.4p below. 

As the Xi's commute, the operator /(x) ■— fixii ■ ■ ■ i Xn) is well defined for any 
polynomial / G K[x]. Then K[x] acts on Span(;B) by 

G K[x] X Span(f?) ^ G Span(B). 

Recall that 1 G B. The mapping 

, . ip: K[x] Span(6) 

^ ^ ^ / - /(x)(i) 

is a homomorphism and, by the following property, 

(2.5) ^{fg) = /(x)(g(x)(l)) = ./(x)(^(5)) V/, g G K[x], 

ker ip is an ideal in K[x]. We now prove that (p coincide on Span(C~'') with the 
projection tt on Span(Z?) along K = ker i/^^ . 

Lemma 2.2. For any element m G , ip{m) — 7r(m). 

Proof. We use induction on the degree of m. If m = 1, we have </3(l) = 7r(l) = 1 
since 1 G .B. Let m 7^ 1 G C"*". As C is connected to 1, to is of the form to — x^toi 
for some toi G C^. By the induction assumption, we have (^(toi) — 7r(TOi). Then, 

ip{m) = (p{ximi) = Xi{^{mi)) = Xj('^("^i)) = XiTi{mi) - k, 

with K E F C K . But we also have 

m ~ Xi mi — Xi{TT{mi) + TOi — 7r(TOi)) = a:j7r(mi) + Xi Ki 

where ki = toi — 7r(mi) G K. We deduce that 

TO — ip{m) + K + Xi Kl = ipijn) + K2 

with K2 = K + Xi Kl G if + n Span(C+). As ki G if and XiKi G Span(C"'"), we deduce 
using (|2.ip that XiKi G K, thus implying K2 G if. As ip{m) G Span(S), it coincides 
with the projection of to on Span(S) along if. □ 

This implies directly: 

(2.6) v{b) = b, p{x,b) = xr{h) "ibeB 'ii = l,...,n, 

(2.7) K{pq) = A(p ^{q)) - A(^(p)^((7)) Vp, q G Span(C+). 
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Lemma 2.3. For all p,q e Span(C+), A{pq) = A{ip{pq)). 
Proof. We first show by induction on the degree of m G that 
(2.8) A(m6) = A{^{mb)) V6 E B. 

The result is obvious if to = 1. Else, as C+ is connected to 1, we can write to = x^toi 
where toi e C+. Using first (|2.7|) and then (j2.6|) . we find: 

A{mb) = A{miXib) = A{miip{xib)) — A{miXi{b)). 

Next, using first the induction assumption and then (|2.5|) . (|2.6[) . we find: 

A(toix^(6)) = A{ip{mlX^{b))) = A(mi(x)(x.,(6))) = A(m(x)(&)) - A(^(m6)), 

thus showing (|2.8p . We can now conclude the proof of the lemma. Let p,q € 
Span(C+). Then, using successively ^J}, (gH), jM]), A{pq) is equal to 

Aip ifiiq)) = A{ip{pip{q))) = A{pix){ip{ip{qm = A{p{x){ip{q))) = Aiip{pq)). □ 

We can now conclude the proof of Theorem ll.51 Let A be the linear operator on 
IK[x] defined by 

A{p) := A{ip{p)) for p e K[x]. 

We show that i?^ is the unique flat extension of i?^^ . 

First, Hj^ is an extension of H^^ since, for all p,q £ Span(C"'"), A{pq) = 
A{ip{pq)) = A{pq) (by Lemma [231). 

Next, we have K = keri?^^ C kcrH^. Indeed, let k E K. By Lemma [2.21 
(P{k) — 7r(K) = 0. Thus for any p E K[x], we have A{pk) = A{ip{pK)) — 
A{p{x){'p{k))) ~ 0, which shows that k E kciH^. 

As is a rewritting family for B and B contains 1, ;B is a generating set of 
K[yi]/{F) and thus dimIK[x]/(i^) < \B\. Set := K[x]/keri?^. Then, as 
F C K C keiHj^, we have dimAjy < dimK[x]/(F) < \B\. On the other hand, 
dim.A^ — TSLiikHj^ > rankTJ^ = ranki/^ = \B\. Therefore, dim^^ — ranki/^ = 

\B\, keiHj^ = (K), Hj^ is a flat extension of H'j^ , and we have the direct sum: 
]K[x] — Span(i3) ©ker Hj^. Moreover, ip{p) is the projection of p S K[x] on Span(i3) 
along kevHj^. Indeed, (p{p) E Span(i3) and p — ip{p) E kei Hj^ for any p E K[x] 
since, for any g G K[x], 

A{pq) = A{^{pq)) = A{p{x){ip{q))), 

A{pip{q)) = A{ip{pip{q))) = A{p{x){ip{q))) = A{pq). 

Finally, if A' E K[x]* is another linear form whose Hankel operator Hfi,i is a flat 
extension of , then ker Hj^ — (K) C ker Ha' ■ This implies that for all p E ]K[x] , 
A'{p) = A'{(p{p)) = A{ip{p)) = A{p). This shows the unicity of the flat extension of 
, which concludes the proof of Theorem 11.51 

3. Applications 

3.1. Application to the flat extension theorem of Curto and Fialkow. 
Theorem 11.21 is in some sense an 'easy' instance of Theorem ll.4l Indeed, under its 
assumptions, one can show existence of a maximum rank principal submatrix of 
Mt-i{y) indexed by a monomial set B connected to 1 which, as noted in the proof 
of Theorem II. 4[ permits to apply Theorem 11.61 
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Proposition 3.1. Let A G (Span(C+ • C+))*, where C := Mn,t-i- //rank H^^ = 
rank H^, then there exists B Q C closed under taking divisions (and thus connected 
to 1) for which rank -ff^^ = rank = \B\. 

Proof. Let M — iA.{ab))g_ f,^c+ denote the matrix of H^* in the canonical bases. 
Consider a total degree monomial ordering ^ of C and let ;B C C index a maximum 
linearly independent set of columns of M which is constructed by the greedy al- 
gorithm using the ordering ^. One can easily verify that B is closed under taking 
divisions (cf. [13]). □ 

The following example shows that, even if C is connected to 1, there may not 
always exist a base B connected to 1 for (which justifies our generalisation of 
Theorem II .61 to kernels of Hankel operators). 

Example. For n = 2, let C = {1, xi, a;iX2} with dC = {x2,xix'^^x\,x\x2} , and let 
A e (Span(C+ • C+))* be defined by A(a;l4) = 1 if j = 0,1, and K{x\x{) = a if 
i — 2,3,4, except K{x\x\) — a^, where a is a scalar with a ^ 1. The associated 
moment matrix has the form 
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and rank i/^ = rank i/^^ =2. As 1 — G ker H'}^ , the only sets indexing a 
column base for are B — {l,a::ia;2} and {xi,xiX2}, thus not connected to 1. 



Combining Theorem 11.51 with Theorem 11.11 we obtain the following extension of 
Corollary O 



Theorem 3.2. Let A G (Span(C"'' ■ C'^))* , where C C A4„ is finite and connected 



to 1. Assume that A is positive and that rank -ff^^ — rank H^. Then the sequence 



y — (A(a))Qgc+.c+ has a representing measure. 

3.2. Application to polynomial optimization. We point out here the relevance 
of the flat extension theorems to polynomial optimization and to the problem of 
computing the real roots to polynomial equations. In this section, we take again 

The truncated moment problem has recently attracted a lot of attention also 
within the optimization community, since it can be used to formulate semidefinite 
programming relaxations to polynomial optimization problems (see |11|). Moreover 
the flat extension theorem of Curto and Fialkow permits to detect optimality of the 
relaxations and to extract global optimizers to the original optimization problem 
(see [8]). Here is a brief sketch; see e.g. [15] and references therein for details. 

Suppose we want to compute the infimum p* of a polynomial p over a semi- 
algebraic set K deflned by the polynomial inequalities gi > 0, . . . , gm > 0. For any 
integer t > deg(p)/2 and such that t > dj :— \deg{gj)/2], consider the program: 

(3.1) p: := inf A{p) s.t. A G (M[x]2t)*, A(l) = 1, A ^ 0, 5, • A ^ (Vj < m). 
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Here, A ^ means that A is positive (i.e., A(p^) > for all p G M[x]t) and the 
localizing conditions 5j • A )^ (i.e. A{gjp^) > for all p G M[x]t_rfj. ) aim to 
restrict the search for a representing measure suported by the set K (cf. [SI lllj). 
Using moment matrices, the program (jS.ip can be formulated as an instance of 
semidefinite programming for which efficient algorithms exist (see e.g. [20l[21]). We 
have: p^ < p* , with equality if H'l^"'* is a flat extension of ".'-t* Jqj. optimum 
solution A to (|3.f p {d := maxjdj). In that case, the atoms of the representing 
measure (which exists by Corollary II. 3p are global minimizers of p over the semi- 
algebraic set K and they can be computed from A [H]. Moreover, they are all the 
global minimizers when iJ^"'* has the maximum possible rank among all optimum 
solutions to the semidefinite program (j3.ip . 

As shown in [12j . the truncated moment problem also yields an algorithmic 
approach to the problem of computing the real roots to polynomial equations gi — 
0, . . . , g„i — (assuming their number is finite). Indeed, this amounts to finding 
all global minimizers to a constant polynomial, say p = 0, over the real variety 
K := {x ^ M" I gj{x) = Vj = 1,...,to}. Consider the semidefinite program 
(|3.ip where the localizing conditions now read gj ■ A ~ \fj. For t large enough, 
the program (j3.ip has a maximum rank solution which is a flat extension and thus, 
as noted above, all points of K can be computed from this solution. See [12] for 
details. 

A concern in this type of approach is the size of the matrices appearing in 
the semidefinite program (|3.ip . In order to improve the practical applicability of 
this approach, it is crucial to derive semidefinite programs involving matrices of 
moderate sizes. For this one may want to consider moment matrices indexed by 
sparse sets of monomials instead of the full degree levels Mn,t- This is where our 
new sparse flat extension theorem may become very useful. It will be used, in 
particular, in [TP] . 

The approach in [12] also permits to find the real radical of the ideal generated by 
the polynomials gi, . . . , gm- Indeed, if A £ (]R[x])* is positive, then the kernel of its 
Hankel operator H^^ is a real radical ideal [M] and, under the conditions of Theorem 
13.21 ker ff^ generates a real radical ideal. These facts explain the relevance of 
moment matrices and Hankel operators to the problem of finding the real radical 
of a polynomial ideal. For instance, this permits to weaken the assumptions in 
Proposition 4.1 of [12] and to strengthen its conclusions; more precisely, we do not 
need to assume the commutativity of the operators Xi 's (as this holds automatically, 
by Lemma [2. ip and we can claim that the returned ideal is real radical (by the above 
argument). 
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